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•  The  discharge  behavior  of  Li-S  cells  in  wide  range  of  rate  constants  for  precipitation  reactions  is  studied. 

•  Dissolution  of  elemental  sulfur  and  precipitation  of  produced  polysulfides  play  an  important  role  in  capacity  performance. 

•  Optimum  sulfur  content  is  obtained  to  maintain  the  maximum  capacity  performance. 

•  Formulating  the  precipitation  reactions  and  electrochemical  reactions  is  needed  to  be  modified  in  future  works. 
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A  sensitivity  analysis  of  a  mathematical  model  of  a  lithium— sulfur  (Li— S)  battery  was  performed, 
focusing  on  the  precipitation  rate  constants  and  sulfur  content,  by  investigating  the  response  of  the 
model  to  the  variation  of  these  parameters  over  a  wide  mathematical  range.  The  necessity  of  rapid 
dissolution  of  elemental  sulfur  and  rapid  precipitation  of  reduced  sulfur  is  discussed  in  detail.  The 
sensitivity  analysis  suggests  modifying  the  reduction  reaction  steps  of  sulfur  can  improve  the  predictions 
of  the  model  of  Li-S  battery.  Furthermore,  an  upper  limit  on  the  sulfur  content  of  the  cathode  exists  to 
ensure  optimal  performance.  Although  the  model  provides  valuable  knowledge  concerning  Li— S  batte¬ 
ries,  a  modification  to  the  assumed  five-step  reduction  of  sulfur  combined  with  the  consideration  of  the 
insulating  nature  of  the  active  material  is  required  to  improve  the  model. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Sulfur  (S)  one  of  the  favorite  active  materials  for  Lithium  (Li) 
batteries  due  to  its  high  theoretical  specific  capacity,  low  cost,  high 
availability,  and  non- toxic  nature  [1].  However,  several  serious 
problems  hinder  the  commercialization  of  rechargeable  Li— S  bat¬ 
teries.  For  example,  the  low  electronic  conductivity  of  sulfur  and 
polysulfides  necessitates  the  addition  of  conductive  materials  to 
the  cathode,  and  the  dissolution  of  polysulfides  leads  to  a  shuttle 
mechanism  that  causes  self-discharging  of  the  cell  [2,3].  Despite 
the  complex  and  not  fully  understood  behavior  of  sulfur  [4-6], 
attempts  have  been  made  to  simulate  Li— S  batteries  to  provide 
insight  into  their  mechanism  of  operation  [7,8]. 

Continuing  the  work  of  Part  I  [9],  this  paper  presents  a  sensi¬ 
tivity  analysis  of  the  model  presented  by  K.  Kumaresan  et  al.  [8  on 
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precipitation  constant  rates.  The  model  considers  a  five-step  elec¬ 
trochemical  reduction  reaction  of  elemental  sulfur  and  also  in¬ 
cludes  the  dissolution  and  precipitation  reactions  of  sulfur  and 
polysulfides  across  the  cell.  However,  the  authors  introduced 
numerous  parameters  whose  values  must  be  determined  through 
the  appropriate  experiments.  In  the  absence  of  such  experiments,  a 
sensitivity  analysis  can  shed  light  on  the  different  aspects  of  the 
model  on  and  clarify  Li— S  cell  behavior.  To  study  the  physical  and 
mathematical  abilities  and  limits  of  the  model,  a  wide  mathemat¬ 
ical  range  is  assumed  for  the  parameters. 

The  behavior  of  the  model  at  various  discharge  current  rates  was 
previously  studied  [9],  and  the  details  of  the  voltage  plateaus  at 
each  current  rate  were  explained.  Furthermore,  a  sensitivity  anal¬ 
ysis  of  the  cathode  conductivity  was  performed,  demonstrating  the 
existence  of  an  essential  threshold  of  conductivity  for  a  working  cell 
[9].  After  discharge,  the  distribution  of  the  precipitants  is  not  uni¬ 
form;  the  degree  of  this  non-uniformity  depends  on  the  discharge 
current  rates  and  cathode  conductivity.  In  addition,  the  precipitants 
diffuse  in  the  separator,  thus  causing  a  fading  cycle  life. 


0378-7753 /$  -  see  front  matter  ©  2014  Elsevier  B.V.  All  rights  reserved. 
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This  part  of  the  sensitivity  analysis  focuses  on  the  rate  constants 
of  precipitation  and  the  sulfur  content. 


2.  Model  development  and  governing  equations 

The  model  is  based  on  the  proposed  model  by  K.  Kumaresan,  Y. 
Mikhaylik,  and  R.E.  White  [8].  Fig.  1  presents  a  schematic  of  the  Li-S 
cell  and  summarizes  the  governing  equations  and  boundary  con¬ 
ditions  of  the  model.  Appendix  B  describes  the  model,  and 
Appendix  A  defines  its  parameters.  The  values  proposed  for  the 
parameters  are  provided  in  Tables  1—4  in  Appendix  C. 

The  reactions  in  the  cell  include  Li  metal  oxidization  at  the 
anode  surface  during  discharge  [8]: 

Li  Li+  +  e-  (1) 

It  is  assumed  that  during  discharge,  the  elemental  sulfur,  which 
is  initially  in  the  solid  phase,  dissolves  in  the  electrolyte  and  then 
goes  through  the  following  electrochemical  reactions  [8]: 


lc  -_1C2- 

2^8(l)+e 

(2) 

?Sl-+e-^2S§- 

(3) 

S26-+e-^S24~ 

(4) 

^Sj-  +e~^±Sl~ 

(5) 

^sj-+e-^±2S2- 

(6) 

The  following  precipitation/dissolution  reactions  are  also 


present: 

^8(s)  ^$8(1) 

(7) 

2Li+  +  Sg  ^Li2Sg(S) 

(3) 

2Li+  +  S4  ^Li2S4(S) 

0) 

2Li+  +  S2  ^Li2S2(S) 

(10) 

2Li+  +  S2  ^Li2S(S) 

HD 

The  details  of  the  model  and  its  governing  equations  can  be 
found  either  in  Appendix  B  or  ref.  [8].  The  current  densities  due  to 
the  electrochemical  reactions  are  given  by  the  Bulter— Volmer 
equation. 


3.  Results  and  discussion 

A  range  of  variations  is  assumed  for  each  parameter. 
This  range  is  not  bounded  by  a  range  of  physical  values.  The 
goal  is  to  determine  the  behavior  of  the  model  system  with 
respect  to  different  situations  and  to  find  a  range  for  the  pa¬ 
rameters  in  which  the  Li-S  cells  are  feasible.  However,  the 
functionality  of  the  model  with  respect  to  these  parameters  was 
found  not  to  be  linear.  Instead,  we  must  choose  a  parameter 
each  time  and  investigate  the  changes  in  the  model  from  vari¬ 
ations  in  this  parameter  while  keeping  the  other  parameters 
constant. 
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Fig.  1.  Schematic  of  a  Li— S  cell  and  a  summary  of  the  governing  equations  and  boundary  conditions.  Repeated  from  Ref.  [9]. 


404 


M.  Ghaznavi,  P.  Chen  /  Journal  of  Power  Sources  257  (2014)  402-411 


The  normalized  current,  due  to  electrochemical  reaction  j  in  the 
cathode  and  defined  by  the  following  expression  9],  is  used  to 
illustrate  the  kinetics  of  the  system: 

x  =  L 

^  =  r~  [  aijdx  02) 

*app  J 

x  =  Ls 

where  J2j= 2  =  1.  In  other  words,  describes  the  overall 

contribution  of  the  electrochemical  reaction  j  in  producing  the 
discharge  current  at  a  given  time.  In  the  ideal  case  of  utilizing  the 
entire  capacity  of  sulfur,  the  discharge  normalized  currents  satisfy 
the  following  relations: 

Tdc 

7^  /  !jdt  =  ^  03a) 

t  =  0 


<13b) 

where  TDC  is  the  total  time  of  an  ideal  complete  discharge  and  Cj  is 
the  portion  of  reaction  j  in  the  total  capacity.  The  responses  of  the 
normalized  currents  with  respect  to  the  applied  current  were  dis¬ 
cussed  in  detail  in  Ref.  [9]. 


3.1.  Rate  constants  for  the  precipitation  reactions 

Each  rate  constant  controls  the  dissolving  and  precipitating 
speed  of  the  related  species.  Although  the  model  predicts  that  the 
precipitation  rate  and  dissolution  rate  behave  symmetrically 
around  the  saturation  concentrations  (see  Equation  (B13)),  far  away 
from  saturation,  they  behave  differently:  there  is  no  limitation  on 
the  precipitation  rates  if  the  concentration  of  the  species  increases, 
but  dissolution  rates  have  an  upper  limit  that  depends  on  the 
related  solubility  product,  rate  constant,  and  solid  volume  fraction. 

This  fact  means  the  precipitated  species  may  not  dissolve  back 
into  the  solution  to  be  involved  in  the  reactions.  In  the  model,  an 
initial  equilibrium  between  the  electrochemical  reactions  and 
precipitation  reactions  is  assumed;  therefore,  changes  in  the  solu¬ 
bility  products  were  avoided  to  conform  to  this  assumption.  Nu¬ 
merical  instability  was  observed  when  this  assumption  was 
ignored. 

The  solubility  of  each  species,  and  thus  their  respective  rate 
constants,  depends  on  the  solvent  used.  Because  the  kinetics  of 
each  reaction  is  influenced  by  the  concentration  of  each  species,  the 
discharge  voltage  will  also  depend  on  the  solvent.  In  this  section, 
the  rate  constant  parameters  are  the  subject  of  study.  Figs.  2  and  3 
present  the  simulation  results  for  0.1  C  and  1  C,  respectively. 

The  behavior  for  the  rate  constant  of  solid  elemental  sulfur,  /<Sg(s) , 
is  simple.  If  the  dissolution  speed  is  fast  compared  to  the  reduction 
rate  of  Ss(i>  a  two-stage  voltage  plateau  exists.  As  the  rate  constant 
decreases,  the  Ss(i)  concentration,  and  thus  the  discharge  voltage, 
becomes  lower.  The  first  discharge  plateau  vanishes,  and  the 


Fig.  2.  Discharge  voltage  plateau  for  different  rate  constants  at  discharge  current  rates  of  0.1  C.  (a)-(e)  In  each  case,  only  the  mentioned  rate  constant  differs  from  the  initial 
assumptions,  (f)  Some  other  cases  for  comparison.  In  each  set,  kk  =  k*k  x  bk,  and  bk  are  equal  to  one  in  “set*”,  {1, 1, 20,  6, 1}  in  “set  1”,  {0.1, 1,  25,  8,  2.5}  in  “set  2”,  and  {0.051, 1,  60, 16, 
250}  in  “set  3”. 
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Fig.  3.  Discharge  voltage  plateau  for  different  rate  constants  at  a  discharge  current  rate  of  1  C.  (a)— (e)  in  each  case,  only  the  mentioned  rate  constant  differs  from  the  initial  as¬ 
sumptions.  (f)  Discharge  voltage  plateau  at  high  C  rates  for  the  case  with  a  high  solubility  of  sulfur. 


contribution  of  the  sulfur  reduction  in  the  discharge  current  per 
time  decreases,  causing  S8(S)  to  last  longer  in  the  system.  As  ex¬ 
pected,  for  equal  ratios  of  /<Sg(s)  //app,  similar  behaviors  for  discharge 
and  retained  discharge  capacity  were  observed  (Figs.  2a  and  3a); 
only  the  sharp  trough  changes  due  to  the  different  ratio  of 

^Li2S(s)/^app. 

For  the  mid-polysulfide  products,  the  behavior  is  more 
complicated,  particularly  on  the  voltage  plateau.  The  effect  of  the 
rate  constant  on  the  capacity  can  be  summarized  as  follows.  To  start 
the  precipitation  of  a  species,  a  rate  constant  threshold  is  required, 
which  depends  on  the  concentration  of  the  species  and  thus  the 
applied  current.  As  rate  constant  exceeds  the  threshold  to  a  greater 
degree,  the  amount  of  the  solid  phase  of  the  species  formed  be¬ 
comes  greater.  As  mentioned  earlier,  rapid  precipitation  does  not 
necessarily  mean  rapid  dissolution,  and  the  precipitated  material  is 
not  involved  in  the  electrochemical  reaction,  thus  leading  to  ca¬ 
pacity  loss.  If  the  rate  constant  is  sufficiently  large  to  reproduce  the 
amount  of  species  dissolved  in  the  electrolyte,  the  retention  ca¬ 
pacity  increases  again.  This  type  of  behavior,  over  a  mathematical 
range,  is  observed  for  /<Li2s8(s)  and  /<s8(s)  (Figs.  2b,  3b  and  3c). 

The  voltage  plateau  behavior  depends  on  more  parameters  and 
is  more  complicated.  At  0.1  C,  increasing  /<Li2s8(s)  by  a  factor  of  400 
causes  the  threshold  of  initiating  precipitation  to  be  exceeded, 
leading  to  some  capacity  loss.  If  the  rate  constant  is  increased 
further,  s|_  will  be  removed  earlier  and  faster,  causing  the  reduc¬ 
tion  of  S8(i)  to  increase  up  to  a  peak;  not  the  peak  of  1^  will  vanish, 
but  the  reverse  oxidization  of  Sg-  will  make  it  negative.  This  phe¬ 
nomenon  occurs  for  /<Li2s8(s)  >  ^u2s8(s)  x  causing  significant 
capacity  loss. 


A  factor  of  104  makes  this  phenomenon  stronger:  as  precipita¬ 
tion  starts,  when  it  passes  its  peak,  1%  becomes  one  and  the  voltage 
increases,  i.e.,  a  new  type  of  deep  trough  forms.  Although  precipi¬ 
tation  has  started,  the  S concentration  remains  constant  while  all 
of  the  elemental  sulfur  reduces.  With  a  decrease  in  the  concen¬ 
tration  of  S8(i),  the  potential  also  decreases,  and  some  of  the  will 
be  reduced.  Only  a  small  portion  of  sulfur  will  reduce  to  S2~,  and 
most  of  the  capacity  will  be  lost. 

^Li2s8(S)  =  ^Li2s8(s)  x  almost  all  of  the  capacity  can  be 
maintained.  Precipitated  materials  dissolve  again,  contributing  in 
electrochemical  reactions. 

A  higher  rate  constant  threshold  is  expected  at  1  C.  Interestingly, 
the  lowest  maintained  capacity  is  much  higher  than  0.1  C,  and  the 
potential  in  the  first  plateau  always  decreases.  The  key  is  the  lower 
dissolution  rate  of  elemental  sulfur  with  respect  to  the  electro¬ 
chemical  reaction  rates.  The  S8(i)  concentration  is  always  low,  and  1% 
cannot  become  one;  therefore,  the  voltage  of  the  first  plateau  re¬ 
duces.  In  contrast,  the  Sg_  concentration  does  not  considerably 
exceed  saturation,  leading  to  slow  precipitation  compared  to  the 
electrochemical  reaction  rate.  Particularly  after  a  period  of  time, 
because  the  Sg_  production  rate  is  low,  its  concentration  remains 
below  saturation,  even  causing  the  precipitated  material  to  re¬ 
dissolve  in  the  electrolyte.  Therefore,  a  cell  discharges  more  even 
for  the  same  ratios  of  /<u2s8(s)/fapp. 

The  same  explanation  applies  to  the  /<Li2s4(s)  parameter  as  well.  If 
this  rate  constant  is  increased,  active  material  is  lost  by  precipita¬ 
tion.  Eventually,  an  increase  in  this  rate  allows  for  re-dissolution, 
thus  retaining  more  capacity.  precipitation  lowers  its  concen¬ 
tration,  causing  a  decrease  in  potential  along  with  a  local  minimum 
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in  1^  and  Jg ,  whereas  1%  presents  a  peak.  The  decrease  in  potential 
allows  for  the  semi-equilibrium  of  reaction  rates. 

At  0.1  C,  only  for  /<y2 S4s  x  500  does  S4-  precipitation  occur 
before  the  second  plateau, leading  to  another  deep  trough  in  the 
potential.  At  this  point,  the  reduction  of  S4-  has  not  yet  started.  As 
the  Sl~  concentration  decreases  due  to  precipitation,  less  over¬ 
potential  is  required  to  reduce  Sg-  and  the  potential  increases, 
leading  to  faster  Sg-  reduction  and  also  to  its  partial  reverse  oxi¬ 
dization  and  the  partial  reverse  oxidation  of  Sg-.  The  increased 
potential  leads  to  a  new  semi-equilibrium  between  the  electro¬ 
chemical  reactions  (J^,  1%  and  1%)  until  the  Sg-  concentration  de¬ 
creases  very  suddenly,  which  causes  a  sudden  decrease  in  potential. 
The  same  behavior  is  observed  at  1  C. 

At  a  rate  constant  increased  by  factors  of  104  and  106,  precipi¬ 
tation  again  occurs  in  the  first  plateau,  causing  another  deep 
trough.  At  factors  of  106  and  higher,  the  discharge  capacity  starts  to 
increase  again. 

Increasing  /<Li2s2(s)  increased  the  amount  of  S2-  precipitation.  As 
the  constant  rate  becomes  larger,  the  precipitation  becomes  faster 
and  the  S2-  concentration  decreases.  Therefore,  the  potential  in¬ 
creases  along  with  an  increase  in  I$j  and  a  decrease  in  l£.  At 
ky2 s2  s  x  the  precipitation  of  S2-  and  S2-  both  begin  simulta¬ 
neously  (see  Fig.  2d).  1%  increases  to  a  very  large  peak,  and  with  an 
increase  in  the  potential,  reverse  oxidation  of  S4-  occurs,  with  1% 
becoming  almost  zero  after  a  sharp  and  relatively  small  peak.  Only 
a  small  portion  of  reduces  to  S2-. 

Increasing  the  rate  constant  by  a  factor  greater  than  50  causes 
the  precipitation  to  start  before  the  second  plateau.  Consequently, 
the  factors  of  105  and  106  cause  a  large  increase  in  the  potential,  and 
no  reduction  of  S2_  occurs.  Therefore,  only  half  of  the  capacity  is 
maintained. 

The  same  behavior  is  observed  at  1  C  (see  Fig.  3d).  For  the  same 
ratios  of  kLi2S  /Iapp,  the  same  discharge  capacity  is  retained.  At¬ 
tempts  to  increase  the  discharge  capacity  by  increasing  the  rate 
constant  failed  because  of  numerical  instability. 

The  precipitation  of  the  last  polysulfide,  S2-,  is  controlled  by 
/<Li2s(s).  Removing  the  last  product  from  the  electrolyte  is  necessary 
to  fully  discharge  the  cell  and  to  avoid  the  cost  of  decreasing  voltage 
due  to  high  product  concentration.  As  long  as  the  rate  constant  is 
above  a  threshold  (which  depends  on  the  applied  current),  the 
product  removal  occurs  perfectly  and  no  change  in  the  shape  or 
behavior  occurs  (see  the  largest  values  of  the  rate  constant  in 
Figs.  2e  and  3e). 

At  0.1  C,  the  only  difference  in  the  discharge  plateau  between 
the  two  highest  values  of  the  rate  constant  is  the  lack  of  the  trough 
(also  in  1  C,  see  Fig.  3e).  At  high  values  of  the  rate  constant,  as  soon 
as  S2-  reduction  begins,  S2-  product  begins  to  precipitate  because 
very  tiny  solid  phase  nuclei  can  initiate  the  precipitation.  However, 
at  ky2s  ,  the  tiny  nuclei  cannot  initiate  the  precipitation,  so  the  S2- 
concentration  exceeds  the  saturation  point  and  the  potential  de¬ 
creases  until  the  very  sudden  start  of  precipitation  removes  the 
product,  at  which  point  the  potential  increases. 

Decreasing  the  rate  constant  by  a  factor  of  0.01  delays  the  pre¬ 
cipitation,  which  causes  the  accumulation  of  S2-  in  the  electrolyte, 
thereby  resulting  in  a  further  decrease  in  the  potential.  The  lower 
potential  causes  enhanced  production  of  as  well,  such  that  it 
starts  to  precipitate;  however,  the  precipitation  does  not  last  long 
because  S2-  suddenly  starts  to  precipitate.  The  removal  of  S2- 
brings  the  potential  back  to  a  higher  level,  along  with  a  very  large 
peak  of  S2_  reduction,  its  reverse  oxidization  to  S4-,  and  even 
reverse  oxidization  to  Sg-.  As  the  rate  constant  continues  to 
decrease,  the  precipitation  becomes  more  delayed  and  these  effects 
become  stronger. 

At  k^  s  x  0.0025,  most  of  the  precipitated  species  are  Li2S2(S); 
therefore,  the  discharge  capacity  only  just  exceeds  half  of  the  total 


capacity.  Note  that  as  the  S4-  concentration  (and  also  Sg-)  de¬ 
creases  due  to  electrochemical  reactions  and  most  of  the  S\~  is 
removed  by  precipitation  and  not  by  the  electrochemical  reactions 
(the  contribution  of  ig  in  total  current  is  weak),  the  potential 
gradually  decreases  to  maintain  the  reduction  of  S4-. 

The  same  behavior  is  observed  at  1  C,  except  that  the  deep 
trough  is  stretched  in  its  time  scale  (Fig.  3e).  This  stretching  occurs 
because  when  the  reduction  of  both  S4-  and  S^-  begins,  elemental 
sulfur  and  the  high  polysulfides  are  still  being  reduced  as  well 
(whereas  at  0.1  C,  no  significant  reduction  of  the  high  polysulfides 
occurs  because  they  are  all  consumed;  see  Fig.  3  in  Ref.  [9]). 
Therefore,  Ig  is  not  as  high  as  in  the  0.1  C  case;  as  a  result,  on  the 
time  scale  of  the  1  C  case,  the  precipitation  is  delayed  longer, 
thereby  leading  to  the  stretched  trough.  When  the  rate  constant  is 
decreased  by  a  factor  of  0.002  or  lower,  precipitation  also  oc¬ 
curs.  At  the  factor  of  0.001,  S2-  precipitation  starts  slightly  before 
S2-  precipitation  and  constitutes  the  major  solid  phase.  Decreasing 
the  rate  constant  by  a  factor  of  10-4  causes  all  of  the  precipitation  to 
be  Li2S2(S)  instead  of  Li2S(S),  resulting  in  a  discharge  of  only  half  of 
the  capacity. 

Three  representative  cases  are  presented  in  Fig.  2f  to  empha¬ 
size  the  nonlinear  behavior  of  the  system.  In  each  set,  kk  =  k*k  x  bk 
and  bk  are  equal  to  one  in  “set*”,  {1, 1,  20,  6, 1}  in  “set  1”,  {0.1, 1,  25, 
8,  2.5}  in  “set  2”,  and  {0.051,  1,  60,  16,  250}  in  “set  3”.  The  only 
difference  between  the  discharge  plateau  of  “set*”  and  “set  1”  is 
the  shorter  second  plateau.  The  higher  rate  constants  for  Li2S4(S) 
and  Li2S2(S)  cause  the  loss  of  capacity.  At  the  end  of  the  discharge, 
the  volume  fractions  are  approximately  0.037  and  0.025  for  £Li2S2(s) 
and  ^Li2s4(S)  *  respectively.  ru2s(s)  drops  from  approximately  0.287  in 
“set*”  to  0.183  in  “set  1”.  In  “set  2”,  decreasing  ks  by  a  factor  of 
0.1  makes  the  first  discharge  plateau  similar  to  that  at  the  1  C 
discharge  rate  (see  Fig.  3).  However,  the  three  last  rate  constants 
lead  to  a  significant  capacity  loss.  Note  that  the  second  plateau 
remains  fairly  flat.  The  volume  fractions  are  approximately  0.039, 
0.037,  and  0.159  for  £Li2s4(s),  £Li2s2(s),  and  £Li2S(s),  respectively.  “Set  3” 
indicates  how  the  rate  constants  allow  only  one  discharge  plateau 
to  occur. 

A  comparison  of  different  discharge  rates  for  the  cell  with  high 
solubility  is  shown  in  Fig.  3f.  Even  for  a  very  high  rate  of  7  C,  all  of 
the  capacity  is  retained.  The  two  plateaus  are  easily  distinguishable. 
In  addition,  the  sharp  curves  become  smooth,  particularly  those  in 
the  deep  trough  that  disappeared.  The  reason  for  this  behavior  is 
the  dependency  of  dissolution  on  the  size  of  the  particle;  for  high  C 
rates,  as  sulfur  particles  become  smaller,  they  dissolve  slower  than 
the  rate  of  the  electrochemical  reactions,  thus  causing  a  smooth 
start  and  end  for  each  electrochemical  reaction  and  thus  a  smooth 
increase  or  decrease  of  the  concentrations  of  the  species.  Conse¬ 
quently,  S2-  starts  to  precipitate  gradually,  not  suddenly,  and  the 
trough  disappears. 

The  location  of  the  precipitated  material  is  very  important.  If  the 
rate  constant  is  high,  the  species  precipitate  inside  the  cathode. 
However,  if  the  rate  constant  is  small,  the  dissolved  species  diffuses 
into  the  separator  and  precipitates  there  as  well.  To  illustrate  this 
effect,  the  volume  fraction  of  Li2S(S)  is  presented  in  Fig.  4,  which 
demonstrates  that  a  further  decrease  in  the  rate  constant  leads  to 
an  increase  in  the  formation  of  the  solid  phase  in  the  separator. 
Such  movement  of  the  active  material  to  the  separator  is  a  cause  of 
capacity  fading  by  cycling.  Considering  the  possibility  of  higher- 
polysulfide  precipitation  makes  this  fading  effect  even  worse.  For 
example,  Fig.  5  shows  the  volume  fraction  of  Li2S2(s)  and  Li2S(S) 
across  the  cell  at  the  end  of  a  discharge.  The  low  rate  constant  of 
/<Li2s2(s)  leads  to  a  uniform  precipitation  across  the  entire  cell. 
However,  a  large  rate  constant  leads  to  the  presence  of  more  solid 
phase  everywhere  compared  to  the  case  of  a  low  rate  constant,  but 
not  distributed  uniformly.  The  volume  fraction  is  higher  in  the 
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Fig.  4.  Volume  fraction  of  Li2S(S),  £u2s(s)  across  the  cell  at  the  end  of  a  discharge  with  a 
discharge  current  rate  of  0.1  C.  The  separator-cathode  interface  is  located  at  x  =  9  pm. 


cathode.  The  dynamics  of  precipitation  was  discussed  earlier  in 
Ref.  [9]. 


3.2.  Sulfur  content 

It  is  desirable  to  increase  the  energy  density  of  a  battery  by 
increasing  the  active  material  content  per  volume  of  the  cathode. 
The  volume  expansion  of  the  sulfur  products  during  the  discharge 
processes  and  the  porosity  loss  are  limitations  that  place  an  upper 
limit  for  the  sulfur  content.  This  section  provides  a  quantitative 
discussion  of  the  cell  behavior  with  respect  to  the  sulfur  content 
based  on  the  discharge  current  rate.  Note  that  based  on  the  model, 
the  entire  cathode-electrolyte  interface  is  considered  to  be  an 
active  surface  for  electrochemical  reactions,  thus  ignoring  the 
insulating  nature  of  sulfur  and  the  precipitated  polysulfides. 
Therefore,  the  negative  phenomenon  of  losing  active  surface  area 
for  electrochemical  reactions  is  not  considered.  For  simplicity,  the 
volume  ratio  of  conductive  material  in  a  cathode  and  its  structure 
(porosity  and  specific  surface  area)  are  assumed  to  remain  constant 
even  if  the  sulfur  content  is  changed. 

Simulations  were  performed  at  the  three  different  discharge 
current  rates  of  0.25,  2.5,  and  25  A  m-2.  Fig.  6  presents  the 
discharge  plateaus  for  various  sulfur  content  levels,  and  Fig.  7 
presents  the  retained  discharge  capacity  as  a  function  of  sulfur 


Capacity  (Ah/g) 

Fig.  6.  Discharge  plateaus  with  different  sulfur  volume  content  at  various  discharge 
current  rates. 

content.  Both  of  the  figures  indicate  that  cells  with  less  than  20% 
sulfur  content  (per  volume  of  cathode)  can  be  fully  discharged, 
even  for  high  discharge  rates  of  approximately  1  C.  At  a  high 
discharge  current  rate  of  25  A  m-2,  the  short  first  plateau  and  a 
vanishing  trough  is  expected,  as  explained  earlier. 

When  the  sulfur  content  is  0.3,  at  higher  discharge  rates,  the 
voltage  plateaus  are  very  similar  to  those  in  Figs.  2d  and  3d  when 
Li2S2(s)  forms,  whereas  at  0.25  A  m“2,  the  discharge  rate  is  similar  to 
that  in  the  plateaus  with  S4-  precipitation  (see  Figs.  2c  and  3c).  At  a 
discharge  rate  of  25  A  itT2,  if  the  sulfur  content  exceeds  30%  of  the 
cathode  volume,  the  first  discharge  plateau  attains  its  full  shape 
because  £S8(s)  is  sufficiently  large  to  keep  the  dissolution  sufficiently 
fast  for  a  relatively  longer  time.  An  almost  constant  Sg(i)  concen¬ 
tration  is  obtained,  and  consequently,  a  sufficiently  high  S con¬ 
centration  in  the  electrolyte  (as  the  reduction  product  of  elemental 
sulfur)  is  achieved,  and  thus,  no  potential  reduction  is  required  to 
involve  S|_  in  the  electrochemical  reaction.  Although  with  the 
sulfur  content  of  0.4,  the  voltage  plateaus  at  high  discharge  rates 
are  similar  to  those  in  previous  simulation  results;  at  low  discharge 
rates,  the  voltage  plateau  has  an  unexpected  shape  and  end.  For  a 
sulfur  content  of  0.5,  all  of  the  discharges  are  observed  to  stop  in  a 
manner  that  is  not  regulated  until  the  end  of  the  electrochemical 
reactions. 

Fig.  7  provides  more  details  of  the  discharge  capacity.  For  a  low 
discharge  rate,  the  discharge  capacity  percentage  decreases  rapidly 


Fig.  5.  Volume  fraction  of  Li2S(S)  (cLi2s(s))  aRd  Li2S2(S)  across  the  cell  at  the  end  of  a 
discharge  with  a  discharge  current  rate  of  0.1  C.  The  separator-cathode  interface  is 
located  at  x  =  9  pm. 
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Fig.  7.  Discharge  capacity  percentage  out  of  the  total  capacity  for  the  cells  with 
different  sulfur  contents  and  at  various  discharge  rates. 
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because  the  reduction  reactions  of  the  different  polysulfides  occur 
more  individually,  i.e.,  a  polysulfide  almost  reduces  completely 
before  the  products  start  the  subsequent  reaction.  Therefore,  the 
concentration  of  high  polysulfides  can  reach  the  saturation  con¬ 
centration  and  begin  to  precipitate.  At  a  discharge  rate  of 
0.25  A  m~2,  even  a  sulfur  content  of  0.225  increases  the  S4-  con¬ 
centration  to  the  saturation  point  such  that  £Li2s4(s)  at  the  cathode 
area  is  almost  0.02  on  average.  This  volume  fraction  increases  to 
approximately  0.334  when  the  sulfur  content  is  0.4,  whereas  rLi2S(s) 
decreases  to  approximately  0.06. 

Increasing  the  sulfur  content  by  over  0.4  causes  a  Li2Ss(S)  solid 
phase  to  form.  This  solid  phase  has  a  sudden  increase  for  a  sulfur 
content  of  over  0.45,  reaching  an  amount  of  approximately  0.27  at  a 
sulfur  content  of  0.5.  Thus,  the  capacity  decreases  at  this  point. 
When  the  sulfur  content  is  0.55,  all  of  the  solid  phase  is  Li2Ss(S>  and 
even  for  0.6,  approximately  0.07  of  the  sulfur  volume  fraction  re¬ 
mains  in  the  cathode  without  being  involved  in  the  reactions.  For 
all  of  these  sulfur  content  levels,  only  an  insignificant  amount  of 
solid  Li2S2(S)  forms  because  of  the  simultaneous  reduction  of  S4- 
and  S2_,  which  prevents  the  S2-  concentration  to  reach  the  satu¬ 
ration  point.  In  the  previous  section,  high  polysulfides  were  shown 
to  precipitate  everywhere  in  the  cell.  Large  amounts  of  precipitated 
material  can  block  the  cathode-separator  interface.  At  a  discharge 
current  rate  of  0.25  A  m-2,  this  blocking  phenomenon  occurs  for 
cells  with  a  sulfur  content  of  0.4  and  greater. 

At  high  discharge  current  rates,  partially  simultaneous  reduc¬ 
tion  of  the  different  polysulfides  prevents  the  saturation  of  one 
specific  high  polysulfide  that  causes  higher  capacity  retention  with 
respect  to  the  low  discharge  currents.  At  a  discharge  rate  of 
2.5  A  nrT2,  the  major  solid  phase  changes  from  Li2S(S)  to  Li2S2(S)  by 
increasing  the  sulfur  content  up  to  0.4  (Li2S(S)  and  Li2S2(S)  become 
approximately  0.162  and  0.354,  respectively).  However,  a  very  small 
amount  of  Li2S4(S)  forms  when  the  sulfur  content  is  0.4  (approxi¬ 
mately  0.02),  but  increasing  the  sulfur  content  to  0.45  increases  its 
volume  fraction  suddenly  (approximately  0.12),  which  causes  an 
even  greater  reduction  in  the  discharge  capacity.  The  details  of  the 
average  volume  fractions  of  the  precipitants  in  the  cathode  area  are 
provided  in  Fig.  8.  At  this  point,  the  cathode-separator  interface 
becomes  blocked  as  well.  The  higher  sulfur  content  leads  to  an 
increase  in  the  Li2S4(S)  precipitates.  The  precipitated  material 
almost  completely  fills  the  pores  of  the  separator  (the  porosity 
becomes  less  than  0.05  in  the  separator  at  a  sulfur  content  of  0.6) 
and  blocks  the  cathode-separator  interface. 
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Fig.  8.  Average  volume  fraction  of  precipitants  across  the  cathode  at  the  end  of  a 
discharge  for  different  discharge  current  rates. 


At  a  high  discharge  rate  of  25  A  m-2,  no  solid  phase  of  Li2S4(S) 
forms,  resulting  in  the  smooth  blue  line  (in  web  version)  in  Fig.  7. 
For  sulfur  contents  greater  than  0.4,  the  cathode-separator  inter¬ 
face  is  blocked.  The  slightly  lower  discharge  capacity  at 
Japp  =  25  A  m”2  compared  to  2.5  A  m-2,  when  the  sulfur  content  is 
between  0.225  and  0.4,  is  related  to  the  relatively  earlier  appear¬ 
ance  of  s\~  at  the  higher  current  rates,  which  causes  slightly  more 
precipitation. 

To  find  the  optimum  sulfur  content  to  gain  the  highest  capacity 
of  the  cell,  the  discharge  capacity  per  surface  area  of  the  cathode  is 
presented  in  Fig.  9.  The  optimum  sulfur  content  is  between  0.225 
and  0.25,  depending  on  the  discharge  current.  However,  consid¬ 
ering  the  movement  of  active  material  to  the  separator  due  to  the 
precipitation  of  Li2Sn>2  precipitations,  one  should  keep  the  sulfur 
content  below  0.225  to  avoid  poor  cyclability. 

The  shapes  of  the  curves  in  the  graphs  shown  in  Figs.  7  and  9 
suggest  a  way  to  distinguish  the  precipitant  as  the  final  product 
of  the  discharge  process.  Each  of  the  points  in  the  figure  which 
exhibits  a  sudden  decrease  in  the  line  that  indicates  the  discharge 
capacity  is  related  to  the  initiation  of  the  precipitation  of  the 
species. 

Note  that  these  results  depend  on  the  coefficients  of  the  pre¬ 
cipitation  rate  constants,  particularly  the  shape  of  the  voltage  pla¬ 
teaus.  However,  the  result  that  at  high  sulfur  contents,  a  higher 
discharge  capacity  is  observed  at  high  discharge  rates,  unlike  at  low 
discharge  rates,  is  a  direct  consequence  of  the  assumed  chain  of 
reduction,  meaning  that  other  sets  of  rate  constants  exhibit  the 
same  results.  At  low  discharge  rates,  the  reduction  of  the  various 
species  does  not  occur  simultaneously,  thus  leading  to  the  satura¬ 
tion  of  high  polysulfides. 

In  carbon-sulfur  composite  cathodes,  as  one  of  the  most  popular 
methods  for  cathode  preparation,  the  melt  sulfur  (or  dissolved  to  a 
solvent)  diffuses  into  the  porous  structure  of  carbon  [3,10].  There¬ 
fore,  theoretically,  the  morphology  of  carbon  structure  does  not 
significantly  change.  The  ratio  of  carbon/sulfur  only  determines  the 
ratio  of  sulfur  volume  to  pore  volume  without  remarkable  change 
in  the  conductivity  of  carbon  matrix.  In  this  case,  the  behavior  of 
the  system  follows  the  results  of  this  section  on  sulfur  content. 
However,  high  amounts  of  sulfur  may  disconnect  carbon  particles 
and  cause  to  decrease  in  the  effective  conductivity  of  carbon  matrix. 
B.  Zhang  et  al.  have  reported  loss  of  specific  capacity  by  increasing 
the  sulfur  content  [10  .  They  used  melted  sulfur  diffusion  into  mi¬ 
cropores  of  carbon.  In  other  methods  of  preparation  for  carbon- 


Fig.  9.  Discharge  capacity  per  unit  surface  area  of  the  cathode  for  different  sulfur 
contents  and  at  various  discharge  rates. 
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sulfur  cathodes,  such  as  mixing  [3]  or  carbon  coating  on  a  surface  of 
sulfur  cathode  [11],  the  morphology  of  the  carbon  matrix  may 
changes  by  different  ratios  of  carbon/sulfur.  However,  through  this 
change  of  morphology,  the  available  surface  area  of  carbon  matrix 
may  not  change  considerably;  the  conductivity  of  carbon  matrix  is 
subjected  to  most  significant  change.  Increasing  carbon  contents 
directly  increases  the  conductivity  of  the  cathode.  We  have  dis¬ 
cussed  the  importance  of  conductivity  of  the  cathode  in  detail  in 
Part  I  [9]. 

4.  Conclusions 

In  this  paper,  the  behavior  of  a  Li-S  cell  mathematical  model 
was  investigated  with  respect  to  a  wide  mathematical  range  of  rate 
constants  for  precipitation  reactions.  The  dissolving  speed  of 
elemental  sulfur  was  observed  to  determine  whether  the  discharge 
voltage  has  either  one  or  two  plateaus:  more  rapid  dissolution  than 
consumption  by  electrochemical  reaction  causes  the  two  flat  re¬ 
gions.  In  fact,  kiJI app  plays  an  important  role  in  the  behavior  of  the 
cell  and  its  capacity.  However,  the  model  indicates  a  strong 
nonlinear  behavior  with  respect  to  this  ratio.  More  specifically, 
there  is  a  “critical  interval”  for  each  rate  constant,  in  which  a  tiny 
variation  in  the  rate  constant  causes  a  large  variation  in  the 
response  of  a  cell,  particularly  in  its  capacity.  The  existence  of  these 
critical  intervals  suggests  that  the  model  requires  modification  in 
its  formulating  of  the  precipitation  reactions. 

Moreover,  the  model  fails  to  reproduce  the  voltage  plateau  of 
the  cases  with  capacity  loss.  In  such  cases,  the  final  products  in  the 
cathode  should  consist  of  Li2S(S),  Li2S2(s>  and  Li2S4(S),  and  some  un¬ 
utilized  sulfur  should  remain  as  well.  Forcing  the  model  to  reflect 
this  situation  reforms  the  voltage  plateau  into  shapes  that  cannot 
be  compared  with  the  experiment  results.  In  other  words,  regard¬ 
less  of  how  well  the  model  works  in  simulating  a  perfect  battery 
that  retains  full  capacity,  it  fails  in  the  simulations  of  typical  bat¬ 
teries  that  lose  capacity.  However,  the  model  still  provides  a 
considerable  amount  of  valuable  information.  For  example,  the 
model  indicates  that  for  slow  precipitation,  more  material  diffuses 
to  the  separator,  leading  to  a  decreased  cycle  life.  Specifically,  this 
decreased  cycle  life  phenomenon  is  stronger  if  the  precipitant  is 
one  of  the  high  polysulfides. 

Most  of  the  phenomena  in  a  cell,  which  are  ignored  in  the  model 
(such  as  the  isolating  nature  of  sulfur  and  polysulfides),  can  reduce 
the  predicted  discharge  capacity;  therefore,  the  model  is  also  used 
to  determine  an  upper  limit  on  the  optimal  sulfur  content.  For  low 
discharge  rates,  the  model  predicts  that  the  discharge  capacity 
percentage  reduces  rapidly  by  increasing  the  sulfur  content 
compared  to  that  at  high  rates.  This  phenomenon  is  related  to  the 
assumed  sequence  of  the  reduction  reactions  chain.  Thus,  some 
modification  to  this  assumption  should  be  implemented  in  future 
models. 

Unfortunately,  all  attempts  to  charge  the  modeled  cell 
failed  due  to  a  low  saturation  concentration  of  Li2S(S)  in  the 
electrolyte.  This  fact  indicates  that  modifications  to  the  model 
are  necessary. 
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Appendices 

A.  Parameters  and  symbols  (primarily  taken  from  ref.  [5]) 

a  specific  surface  area  of  the  cathode 

a0  initial  value  of  a 

b  Bruggeman  coefficient 

q  concentration  of  species  i(i  =  Li+ ,  S8(1) ,  S8_ ,  Sg- ,  S4- ,  , 

S2~,  and  A~  (anion  of  the  lithium  salt  used  in  the 
electrolyte)),  mol  nrr3 

Q.ref  reference  concentration  of  species  i,  mol  m-3 

Cj  portion  of  reaction  j  in  the  total  capacity  for  an  ideal 

complete  discharge 

Dit 0  diffusion  coefficient  of  species  i  in  the  bulk  medium, 
m2  s-1 

Di  diffusion  coefficient  of  species  i  in  the  porous  medium, 

2  -i 

m  s 

F  Faraday  constant,  C  equi-1 

ij  current  density  due  to  reaction  j,  A  m-2 

ij, ref  exchange  current  density  of  the  electrochemical  reaction  j 

at  the  reference  concentrations,  A  m-2 
ii  superficial  current  density  in  the  liquid  phase,  A  m-2 

is  superficial  current  density  in  the  solid  phase,  A  m  2 

/app  applied  current  density,  A  m-2 
IF  normalized  current  due  to  electrochemical  reaction  j 

Ksp,i<  solubility  product  of  precipitate  k 
l<k  rate  constant  of  precipitate  l< 

L  thickness  of  the  cell,  m 

Ls  thickness  of  the  separator,  m 

Ni  superficial  flux  of  species  i,  mol  m2  s-1 
rij  number  of  electrons  transferred  in  electrochemical 

reaction  j 

Pij  anodic  reaction  order  of  species  i  in  electrochemical 
reaction  j 

qij  anodic  reaction  order  of  species  i  in  electrochemical 
reaction  j 

R  gas  constant,  J  mol-1  K-1 

Ri  production  rate  of  species  i  due  to  precipitation  reactions, 

mol  m3  s-1 

Rk  rate  of  precipitation  of  solid  species  /<,  mol  m3  s”1 
rz  production  rate  of  species  i  due  to  electrochemical 

reactions,  mol  m3  s-1 

si  j  stoichiometric  coefficient  of  species  i  in  electrochemical 

reaction  j 

T  temperature,  K 

TDc  Total  time  of  an  ideal  complete  discharge,  s 
t  time,  s 

Uj  Standard  Open  Circuit  Potential  (OCP)  of  electrochemical 
reaction  j 

Ujt ref  OCP  of  electrochemical  reaction  j  at  reference 
concentrations,  V 

Vk  molar  volume  of  the  precipitate  k,  m3  mol-1 
Zi  charge  number  of  species  i 

aaj  anodic  transfer  coefficient  of  reaction  j 
aCj  cathodic  transfer  coefficient  of  reaction  j 

e  porosity  of  the  separator  and  cathode 

£k  volume  fraction  of  precipitate  k  in  the  separator  and 

cathode 

cpi  Potential  in  the  liquid  phase,  V 

cps  Potential  in  the  solid  phase,  V 

jitk  number  of  ionic  species  i  produced  by  dissociation  of 
precipitate  k 

Pj  overpotential  for  electrochemical  reaction  j 
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a  effective  conductivity  of  the  solid  phase  of  the  cathode, 

Ni  =  0 

?  morphology  parameter 


Because  the  charge  can  enter  or  leave  the  liquid  phase  only  by 
electrochemical  reactions,  the  following  equations  apply  at  the 
liquid/solid  interface: 


B.  Governing  equations 

The  model  includes  the  electrochemical  reactions  given  by 
Equations  (1)— (6)  and  the  precipitation/dissolution  reactions  given 
by  Equations  (7) — (11 ).  In  a  porous  medium,  the  governing  equation 
for  the  material  balance  of  an  individual  species  is 

r )rr . 

_1=-V.  JV.+r.-R;  (Bl) 


VI,  =  ajfij  (BIO) 

j 

and 

V-is  +  V-i;  =  0  (Bll) 

The  rate  of  consumption  or  production  of  species  i  due  to  pre¬ 
cipitation/dissolution  is  related  to  the  rate  of  precipitation/disso¬ 
lution  reaction  k  by 


where  the  flux  is  given  by 

Y  =  -DiVCi-Zi^FCjVq,,  (B2) 

The  diffusion  coefficient  for  species  i,  Du  is  corrected  based  on 
Bruggeman’s  expression  for  porosity  and  tortuosity:  D\  =  DIt oeb.  The 
rate  of  production/consumption  of  species  i  due  to  the  electro¬ 
chemical  reactions  can  be  written  in  the  form 


Ri  =  (B12) 

k 

where  the  rate  of  precipitation  of  solid  species  k  is  assumed  to  be 
governed  by  the  following  kinetic  equation: 


(B13) 


ri 


^ n,F 


(B3) 


The  volume  fraction  of  the  precipitate,  k,  is  a  function  of  time 
and  is  given  by 


where  the  stoichiometric  coefficients  are  given  in  Table  1,  Appendix 
C.  The  specific  surface  area  of  the  cathode  varies  due  to  the  pre¬ 
cipitation/dissolution  of  the  various  lithium  sulfide  species  and  is 
assumed  to  be  governed  by  the  empirical  expression: 


a  =  a0 


^initial 


(B4) 


where  the  empirical  parameter,  f,  is  assigned  a  value  of  1.5.  The 
Butler— Volmer  equation  yields  the  current  density  due  to  the 
electrochemical  reactions: 


ij  —  *0jref 


{ 


x  exp  - 


(B5) 


where  the  overpotential  for  the  reaction  j  is 

Vj  =  <Ps~  <Pl~  Uj, ref  (B6) 

The  terms  =  Sy  refer  to  anodic  species  and  qij  =  -Sy  refer 
to  cathodic  species.  The  open-circuit  potential  for  reaction  j  is 
given  by 


Ui 


’j,  ref 


=  uf- 


RT 

njF  S  Sl’ 

J  i 


In 


Q,ref 

Tooo 


The  liquid  phase  current  density  is  given  by 


i,  =  Fj>N, 

i 

The  solid  phase  current  density  follows  Ohm’s  law 


(B7) 


(B8) 


d£k  V  R' 

~dt  -  VkRk 

(B14) 

Therefore,  the  porosity  variation  with  time  is 

k 

(B15) 

The  boundary  conditions  of  the  model  are  as  follows.  At  the 
interface  of  the  cathode  and  the  current  collector,  x  =  L,  the  BCs  are 

Ni  =  0 

(B16a) 

h  —  Japp 

(B16b) 

ii  =  0 

(B16c) 

At  the  cathode-separator  interface,  x  =  Ls, 

seperator  ~  ^i, cathode 

(B17a) 

U,  seperator  =  k, cathode  =  Japp 

(B17b) 

is  =  0 

(B17c) 

At  the  surface  of  anode,  x  =  0, 

o 

II 

s? 

(B18a) 

Nf  =  0 

(Bl 8b) 

Ni  =  U/F 

(B18c) 

ii  =  FN'l 

(B18d) 

(B9) 


The  governing  equations  can  be  solved  numerically. 


is  =  -aV(ps 
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C.  Parameter  tables  for  reference  values 


Table  1 

Stoichiometric  coefficients,  Sy  [8]. 


Sij 

Reactions  given  by  equation  (j) 

Species  (i) 

1 

2 

3 

4 

5 

6 

Li+ 

-1 

0 

0 

0 

0 

0 

S8(l) 

0 

-1/2 

0 

0 

0 

0 

si- 

0 

1/2 

-3/2 

0 

0 

0 

q2- 

->6 

0 

0 

2 

-1 

0 

0 

c2- 

0 

0 

0 

3/2 

-1/2 

0 

c2- 

b2 

0 

0 

0 

0 

1 

-1/2 

s2- 

0 

0 

0 

0 

0 

1 

A~ 

0 

0 

0 

0 

0 

0 

Table  2 

Kinetic  and  thermodynamic  properties  [8]. 

Reaction  (j) 

WA  m“2) 

<j 

nj 

uf 

1 

0.394 

0.5 

0.5 

1 

0.0 

2 

1.9719 

0.5 

0.5 

1 

2.39 

3 

0.019719 

0.5 

0.5 

1 

2.37 

4 

0.019719 

0.5 

0.5 

1 

2.24 

5 

1.972  x  10"4 

0.5 

0.5 

1 

2.04 

6 

1.972  x  10"7 

0.5 

0.5 

1 

2.01 

Table  3 

Transport  properties  and  reference  concentrations  [8]. 

Species  (i) 

Zi 

D?0(m2  s-1) 

Ci\ref  (mo1  m~3) 

Li+ 

+1 

1  x  10”10 

1001.08 

S8(l) 

0 

10  x  10"10 

19.0 

c2- 

->8 

-2 

6  x  lO"10 

0.1832 

c2— 

->6 

-2 

6  x  lO"10 

0.3351 

c2— 

->4 

-2 

X 

o 

o 

0.02146 

-2 

1  X  10"10 

5.999  x  10  7 

S2- 

-2 

1  x  10“10 

9.94  x  lO"10 

A~ 

-1 

4  x  10"10 

1000 

Table  4 

Separator  and  cathode  parameters  [8]. 


Parameter 

Separator 

Cathode 

Thickness  (m) 

9  x  10"6 

41  xlO-6 

^initial 

0.37 

0.778 

£S8(s)  Initial 

1  x  10-12 

0.160 

£Li2S8(S)  Initial 

1  x  10"6 

1  x  10"6 

£Li2S4(S)  Initial 

1  x  10"6 

1  x  10"6 

£Li2S2(s),  initial 

1  x  10"6 

1  x  10"6 

£Li2S(S) , initial 

1  x  10-7 

1  x  10"7 

a0 

- 

132,762 

Table  5 

Parameters  for  the  precipitation  reactions,  either  assumed  or  from  Ref.  [8]. 


Precipitate  Rate  constant  (/<£) 

(k) 

Solubility  product  (K,* ) 

Molar  volume 

(Vk(m3mo\~')) 

S8(s) 

1.0  s  1 

19.0  mol  m  3 

1.239  x  10"4 

Li2S8(s) 

1  x  10"11  m6  mol2  S_1 

183400  mol3  m~3 

1.361  x  10"4 

Li2S4(s) 

9.98  x  10-12  m6  mol2  S_1 

21480  mol3  m-3 

7.415  x  10"5 

Li2S2(S) 

9.98  x  10"9  m6  mol2  S"1 

0.6006  mol3  m-3 

4.317  x  10"5 

Li2S(S) 

6.875  x  10-5  m6  mol2  S_1 

1  9.95  x  10-4  mol3  m  3 

2.768  x  10"5 
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